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EXECUTIVE  SUMMARY 


This  Is  a  report  of  a  study  conducted  under  the  Port  Systems  Project  for 
the  Naval  Facilities  Engineering  Command.  An  objective  of  the  Port  Systems 
Project  is  to  investigate  methods  of  predicting  the  utility  requirements  of  the 
port  facilities  that  must  be  met  to  support  future  f  leets.^Theanalysis  and 
results  presented  in  this  report  are  a  refinement  and  extension  of  the  work 
reported  in  A  Methodology  for  Statistical  Analysis  of  Port  System  Requirements. 
The  fundamental  advance  offered  by  this  study  is  that  predictions  can  be  based 
on  recorded  data  rather  than  DM-25  capacity  information.  _ - — 


- -  This  report  starts  with  a  brief  overview  of  the  first  phase  of  the  method¬ 
ology  development  and  continues  with  the  technical  approach  and  the  principles 
of  statistical  theory  applied.  Examples  are  used  to  relate  these  principles  to 
the  analysis  of  electric  power  requirements  as  an  aid  in  understanding  the 
analytical  process  Inherent  In  the  Methodology.  An  overview  of  the  analytical 
procedures  that  constitute  the  basis  of  the  Port  Systems  Requirements  Predic¬ 
tion  Methodology  is  presented.  Computer  programs  used  to  calculate  statistics 
are  appended.  Finally  the  results  achieved  are  presented. 

The  unique  feature  of  this  approach  to  port  utility  prediction  is  that  it 
identifies  and  accommodates  the  contribution  made  by  correlating  the  time  usage 
of  electric  power  consumption  data  among  different  ship  classes.  This  results 
in  non-zero  correlation  coefficients,  which  have  an  affect  on  the  quality  of 
the  prediction.  It  is  shown  that  as  the  number  of  ships  using  electrical  power 


increases,  the  accuracy  of  the  prediction  is  reduced. 
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I.  INTRODUCTION 


1.0  Background* 

The  United  States  Navy  is  presently  undergoing  many  changes.  The  fleet  is 
being  modernised  to  include  600  active  ships  equipped, with  the  most  modern 
systems  available.  The  future  port  requirements  of  these  ships  will  be 
different  from  those  of  today's  mix  of  ships.  Port  Master  Plans  under 
development  recognise  that  port  facilities  and  services  must  be  upgraded  to 
meet  the  needs  imposed  by  a  .changing  fleer  with  changing  operating  procedures. 
It  is  reasonable  to  conduct  research  into  existing  planning  procedures  in  order 
to  Identify  modern  methods  which  can  be  used  to  reflect  fleet  needs  into  the 
port  master  planning  process.  Current  procedures  are  cumbersome  and  allow  very 
little  interaction  for  conducting  tradeoff  studies.  This  study,  focused  on 
analyzing  and  predicting  electric  power  requirements,  is  a  first  step  in 
developing  these  analytical  procedures.  It  was  conducted  under  the  auspices  of 
the  Port  Systems  Project  for  the  Naval  facilities  Engineering  Command. 


1.1  Objective  of  the  Port  Systems  Requirements  Prediction  Methodology  (PSRPM). 

The  objective  of  PSRPM  is  the  development  of  a  statistically  based 
procedure  for  analyzing  and  predicting  port  utility  requirements.  Eventually, 
it  is  envisioned  that  these  procedures  can  be  used  as  the  basis  of  an  inter¬ 
active  computer  program  which  would  permit  an  efficient  analysis  of  alternative 
means  of  upgrading  port  utility  systems* 

This  study  was  conducted  in  two  pares.  The  described  in  Summary 
Report:  An  Application  of  the_  Port  Systems . Requirements  Prediction  Methodology 
(PSRPM)  Using  Electrical  Measurements  at  Sewel'  'a  Point  (12  -  76  to  1  -  77), 
established  the  feasibility  and  the  basic  methodology  of  a  statistically  based, 
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computet  aided  forecasting  process.  F.>r  convenience,  in  this  report,  the  first 
portion  of  this  scud?  will  be  referred  to  as  Phase  I.  The  extension  of  this 
study  is  the  subject  of  this  report  and  it  vlli  be  referred  to  as  Phase  II. 

1*2  Purpose. 

Th  purpose  of  this  report  is  to  document  the  work  accomplished  to  date  in 
the  developaant  of  the  PSRPM  Phase  II. 

1.3  Scope. 

This  document  includes  a  summary  of  the  work  previously  accomplished  in 
PSRPM  (Phase  1)  as  veil  as  a  detailed  description  of  the  work  accomplished  in 
Phase  II.  A  description  of  the  concepts  supporting  the  methodology  and  a 
review  i  the  statistical  procedures  used  in  the  methodology  are  presented 
first.  Next,  the  application  of  these  statistics  and  the  computer  programs 
associated  with  the  statisticcl  application  are  included.  Finally,  conclusions 
and  recommendations  based  on  the  results  of  the  Phase  II  investigation  are 
discussed. 
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2.0  Objective 


fi 

$ 

A 

i 

I 

* 

i... 

'V. 

V  s' 

S 


I 


i 

I 


f 

i 

t 


* 

£ 


a 

*3 


1  W  ■ 


11.  REVIEW  OF  PHASE  I 


The  objective  of  Phase  1  waa  to  develop  tha  foundation  of  a  statistically 
baaad  methodology  which  could  be  uaad  aa  a  tool  by  port  designers  to  predict 
port  utility  requirements. 

2.1  Technical  Approach  to  Phase  I. 

Phaae  1  aaaumed  that  the  requirement  for  a  port  supplied  utility  to  a 
group  of  ships  could  be  estimated  by  statistical  methods.  If  the  utility 
requirement  for  each  class  was  known  and  the  probability  that  the  ship  would 
need  this  utility  service  was  given,  then  a  statisical  prediction  of  the  total 
utility  requirement  can  be  made.  Phaae  I  used  demand  data  from  NAVFAC  DM-25 
and  ship  populations  and  cold  iron  probabilities  were  derived  from  OPNAV  INST 
3111.14U  and  NAVSEC  Report  6139-72-2  (1977),  respectively.  This  information 
was  used  as  input  to  create  a  computer  simulated  port  activity  model.  This 
model  was  then  used  to  compute  the  probability  that  a  selected  ship  mix  would 
create  a  demand  for  electric  powev  greeter  than  a  port's  deliverable  capacity. 
This  probability  is  plotted  against  capacity  to  yield  a  goal  product  on  the 
form  of  Figure  1.  Figure  2  models  the  Phase  I  activity  flow  which  leads  to  the 
construction  of  the  graphical  representation  that  shows  the  probability  of 
failure  to  meet  the  supplied  demand. 

Electric  power  was  chosen  as  the  initial  utility  to  be  studied  because  it 
can  be  easily  measured  and  the  data  can  be  easily  formatted  for  use  in  a 
mathematical  model. 
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Figure  2.  PSRPM  Phase  I  Activity  Flow 

2.2  Assumptions  -  Phase  I. 

The  first  phase  of  PSRPM  development  applied  the  following  simplifying 
assumptions: 

•  The  only  utility  analyzed  was  electric  power,  and  each  ship  was  assumed 
to  be  drawing  its  full  design  capacity  for  shore  power  in  accordance 
with  DM-25 . 

•  No  consideration  was  given  to  the  time  of  day  variations  that  would  be 
expected  for  electric  power  demand. 

•  Ship  class  mix  was  determined  by  Monte  Carlo  (random  number  generation) 
simulation. 

•  The  number  of  ships  at  cold  iron  in  each  ship  class  is  assumed  to  be 
binomially  distributed  and  independent  of  other  classes. 

•  Resource  demands  per  ship  do  not  correlate  with  the  number  of  ships 
present. 

The  logic  associated  with  the  computational  process  associated  with  PSRPM  is 
presented  in  Figure  3. 
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2.3  Output  ol  Phase  I. 

The  primary  output  of  Phase  I  is  a  mathodology  to  produce  a  graphical 
representation  of  the  probability  that  the  demand  for  electric  power  will 
exceed  a  given  capacity.  The  general  format  of  this  product  was  presented  in 
Figure  1.  The  difference  between  the  curves  in  the  Y  direction,  with  X 
constant,  shows  the  change  in  the  probability  of  failure  to  supply  a  required 
demand  that  will  occur  if  the  resource  is  not  upgraded  for  future  demands.  The 
difference  in  the  X  direction,  with  Y  constant,  shows  to  what  level  of  output 
the  resource  must  be  upgraded  to  maintain  the  current  (if  satisfactory;  failure 
of  probability  to  meet  the  total  demand  level. 

Figure  4  is  an  example  of  the  actual  Phase  I  product. 

2.4  Value  of  Phase  I. 

Phase  I  of  the  PSRPM  demonstrated  that  the  capability  exists  to  take 
information  regarding  a  port's  utility  demand  and  produce  a  tool  usable  in  the 
forecasting  of  future  utility  demands.  The  study  utilized  design  capacity 
information  and  applied  proven  statistical  techniques  to  develop  the  methodo¬ 
logy.  The  power  of  the  PSRPM  lies  in  the  application  of  its  fundamental 
product  i.e.,  P  (Demand  > Capacity)  versus  Capacity  curves.  Specifically,  Navy 
planners  and  port  system  designers  can  assemble  alternative  ship  mixes, 
representing  future  fleet  scenarios,  and  use  the  PSRPM  to  generate  its  goal 
product  for  each  alternative.  These  products  can  then  be  used  to  predict  and 
compare  the  power  demand  of  each  alternative  ship  mix.  This  is  a  major 
advantage  over  the  DM-25  method  of  establishing  pier  capacity  because  it  can 
represent  ship  types  not  yet  constructed. 
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HI.  PSRPM  -  PHASE  II 

3.0  Introduction 

Phase  II  of  the  PSRPM  provides  a  refinement  to  the  methodology  by  using 
actual  recorded  electrical  demand  data  as  input  to  the  computational  process 
in  place  of  DM-25  Present  Demand  data  (design  electrical  capacity).  Because 
there  is  variation  in  the  demand  when  surveyed  data  is  used,  the  statistical 
variance  of  each  ship  mix  becomes  increasingly  important.  Previous  studies 
have  based  their  prediction  of  overall  power  demand  for  an  ensemble  of  ships  on 
the  calculation  of  individual  ship  means  and  individual  variances.  PSRPM  Phase 
II  estimates  the  overall  power  demand  of  a  group  of  ships  by  calculating  the 
individual  ship  means  and  the  total  ship  ensemble  variance.  The  governing 
basis  for  the  ability  to  predict  a  ship  mix's  utility  is  the  Central  Limit 
Theorem. 

3.1  Objective 

Phase  II  of  the  PSRPM  had  the  objective  of  applying  the  modeling 
techniques,  developed  in  Phase  I,  to  the  analysis  of  electrical  demand  data 
obtained  from  Sewell's  Point  in  1977. 

3.2  Scope  of  Phase  II 

Phase  II  of  the  PSRPM  was  written  so  the  modeling  techniques  used  in  Phase 
I  could  be  applied,  with  modification,  to  analyze  actual  electrical  power 
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utilization.  The  principle  differences  between  the  two  phases  are  listed 
below: 

*  Phase  II  analyzed  the  actual  three  phase  electrical  power  used  by 
ships  in  port  at  cold  iron.  Phase  I  assumed  each  ship  drew  its 
full  design  capacity  for  shore  power  in  accordance  with  DM-25. 

*  Phase  I  did  not  take  into  consideration  any  time  of  day  variations 
which  occur  in  real  life.  The  Phase  II  analysis  included  a  time 
variance  factor  based  on  the  hourly  electrical  measurements  to 
account  for  day  to  day  fluctuations  in  power  consumption. 

*  Phase  I  assumed  that  the  resource  demands  per  ship  were  not  related 
to  the  number  of  ships  present.  Phase  II  recognizes  that  the 
electrical  demand  of  one  class  of  ships  can  in  fact  be  correlated 
with  demand  for  other  classes.  This  affect  is  determined  by  the 
calculation  of  interclass  and  intraclass  covariability  of 
electrical  power  demand  as  a  function  of  the  total  load. 

*  The  Phase  II  analysis  is  based  on  one  day  "snapshots"  since  the  mix 
of  ships  in  the  port  varied  from  day  to  day. 

The  output  of  Phase  II  is  identical  in  format  to  that  of  Phase  I. 
However,  the  affect  of  using  measurement  data  will  be  seen  to  cause  an  overall 
reduction  in  the  predicted  power  demand  for  the  case  analyzed.  Figure  5 
displays  the  activity  flow  in  Phase  II. 
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Figure  5.  PSRPM  Phase  II  Activity  Flow 


3.3  Phase  II  Process 

The  technical  objective  of  the  Phase  II  computational  process  was  to 
derive  the  distribution  of  electic  power  demand  for  a  group  of  ships  of  various 
classes.  The  process  for  obtaining  this  objective  is  outlined  in  Figure  6. 
The  objective  class  distributions  have  been  assumed  to  be  normal  and  can  be 
defined  by  knowing  a  mean  KVA  utilization  and  the  variance  in  the  utilization. 
Data  for  computing  the  distributions  was  obtained  at  Sewell's  Point.  From  this 
data,  ship  specific  mean  hourly,  mean  daily,  and  hourly  variance  in  KVA  usage 
was  determined.  Next,  the  covariance  associated  with  every  combination  of 
ships  was  calculated.  This  information  constitutes  what  is  necessary  to 
generate  a  port  mix  of  ships  or  realization. 

A  realization  is  generated  by  knowing  the  probability  of  a  ship  class 
being  at  cold  iron  and  the  number  of  ships  in  each  class  asigned  to  the  port. 
A  Monte  Carlo  Simulation  determines  the  ship  mix  or  ensemble  of  ship  classes 
which  make  up  each  of  the  100  realizations.  Normal  probability  distributions 
representing  the  KVA  demand  per  class  were  constructed  based  on  the  class  me  in 
and  the  ensemble  variance.  Ten  different  demand  levels  were  determined  for 
each  class  using  a  Monte  Carlo  Simulation  and  the  probability  distributions. 


The  resulting  1000  port  realizations  of  demand  were  used  to  construct  a  cumula¬ 
tive  probability  curve  displaying  the  Probability  of  Demand  Exceeding  Capacity 
versus  Capacity.  This  P(D>C)  vs  C  curve  and  its  Interpretation  fora  the  basis 
of  the  PSRPM, 

3.4  Assumptions 

Enhancements  to  the  second  phase  of  the  PSRPM  analysis  are  based  on 
changes  aade  to  the  initial  Phase  I  assumptions.  The  new  assumptions  are: 

°  Ship  class  electrical  demand  is  based  on  actual  measurements  rather 
than  DM- 25  values, 

°  The  hourly  variations  In  electrical  power  usage  are  examined,  and, 

°  The  resource  demands  per  ship  are  correlated  with  the  number  of  ships 
present. 


3.5  Sewell *s  Point  Data 

The  ship  class  electrical  demand  information  which  was  used  in  Phase  II 
was  obtained  from  measurements  taken  at  Sewell' i  Point,  Norfolk,  Virginia  in 
December  1976  through  January  1977.  Hourly  samples  of  amperes  were  recorded  at 
450  -  470  volts  AC  for  each  ship  present  during  the  measurement  time  frame. 
The  amperes  were  converted  to  kilovolt  amperes  (KVA)  because  components  of 
electrical  power  transmission  ate  usually  rated  in  KVAs. 


Figures  7a  and  7b  are  time/date  tables  of  the  electrical  measurements  for  ACR-1 
and  CGN-38  ships.  This  data  is  representative  of  the  data  taken  for  the  ships 
present  each  day  at  Sewell's  Point. 
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IV.  PHASE  II  -  TECHNICAL  BACKGROUND 
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4.0  Introduction 

The  PSRPM  was  developed  aa  a  prediction  tachnlqua  for  uaa  by  port  planners 
In  tha  updating  of  praaant  porta  and  in  tha  conatructlon  of  new  onaa.  Tha 
baais  of  thla  prediction  methodology  ia  tha  determination  of  the  probability 
diatribution  function  for  the  electrical  demand  of  a  group  of  ahlpa.  Thia 
distribution  can  then  be  used  to  calculate  tha  probability  of  the  demand 
r , receding  the  facility's  capacity.  Thla  is  the  basis  for  predicting  the  power 
demand  of  a  given  group  of  ships. 

A  distribution  function  is  defined  by  its  mean,  which  locates  the  center, 
and  by  its  variance  which  determines  the  shape.  The  variances  and  means  of 
electrical  demand  of  several  ship  classes  were  calculated  from  the  data 
recorded  at  Sewell's  Point  1976.  Since  the  variance  describes  tha  shape  of  the 
distribution,  its  accuracy  determines  the  quality  of  the  prediction  which  can 
be  made  from  the  PSRPM.  This  study  shows  that  there  is  a  correlation 
associated  with  the  ships'  time  of  electrical  usage.  The  existence  of  a  non 
zero  correlation  coefficient,  by  definition,  establishes  a  covariance  term 
which  must  be  included  when  calculating  the  variance  in  electrical  demand 
caused  by  a  group  of  ships.  The  contribution  of  this  term  to  total  demand  is 
not  intuitively  obvious.  It  will  be  seen  that  the  covariance  term  will  Impact 
the  accuracy  of  a  prediction  system  based  on  statistical  processes. 

The  methodology  developed  herein  utilizes  particular  statistical  techni¬ 
ques  which  provide  significantly  more  information  than  simply  the  maximum 
electrical  capacity  required  at  Naval  ports.  This  section  is  devoted  to  the 
explanation  of  each  of  these  techniques  and  their  relevance  to  the  PSRPM. 
Included  are  discussions  on  the  mean,  the  variance,  the  covariance,  the 
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standard  deviation,  and  the  correlation  coefficient.  Also  included  is  a  review 


of  the  conservative  factor  (K)  which  was  developed  as  part  of  Phase  II  of  the 
PSRPM.  The  details  of  this  factor  are  addressed  in  A  Comparison  of  Two  Methods 
for  Estimating  the  Electrical  Supply  Capacity  Required  for  Piers  and  Drydocks. 


i 


4. 1  The  Normal  Distribution 

Figure  5  displays  a  graph  of  a  normal  density  function,  sometimes  called 
the  Standard  Normal  Curve.  The  shape  of  the  normal  curve  is  defined  by  the 
magnitude  of  the  mean  and  the  standard  deviation,  or  the  square  root  of  the 
variance.  Plus  and  minus  one  standard  deviations  from  the  mean  ifines  the 
inflection  points  of  the  curve.  Within  this  area  lies  68.27%  of  the 
observations. 


8 


a 


Figure  8.  Standard  Normal  Curve 

When  the  variance  in  a  set  of  data  changes,  the  overall  shape  of  the  curve 
changes.  Figure  9  shows  the  effect  of  increasing  variance  on  the  shape  of  the 
distribution  curve.  Phase  II  is  interested  in  the  shape  of  the  normal  curve 
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when  the  variance  increases.  The  change  in  the  shape  of  the  curve, 
particularly  in  the  wings,  reveals  the  degree  of  prediction  accuracy  inherent 
to  the  methodolgy.  Also  the  new  shape  redefines  the  maximum  electrical  demand 
for  a  given  mix  of  ships. 

4.2  Mean 

One  parameter  of  an  entire  group  or  population  that  is  of  interest  is  the 
population  mean,  p.  ,  because  it  is  one  method  used  to  describe  the  population's 
center.  If  every  item,  or  ship  in  the  case  of  this  report,  were  totally 
accessible  and  measured  its  mean  would  be  computed  by  the  formula: 


where 


N 


(1) 


T  Jjk 

i-1 


fl  m  population  mean 


N 

& 

N 


sum  of  all  values  of  the  variable  of  interest  for  the  whole 
population 

Number  of  elements  in  the  population 


*■  * 


r-* 
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Statistical  inference  about  a  population  mean  requires  information 
obtained  from  samples  of  the  population.  In  the  case  of  the  Sewell's  Point  data 
the  information  gathered  are  samples  of  the  population.  The  sample  average,  "x, 
which  is  analogous  to  a  population  mean  is  defined  as: 


x 


n  X- 

E— 

i=l  n 


where  x  «  value  of  the  average 

n 

£x^  ■  sum  of  samples 

i=l 

n  ■  number  of  samples 


(2) 


4,2.1  Summation  of  Means 

If  two  or  more  3ets  of  data  of  equal  sample  size  have  been  collected,  the 
mean  of  the  sum  (xj  +  x2)  is  equal  to  the  sum  of  the  means.  Example  1  is  an 
example  of  the  application  of  this  statement.  Data  was  extracted  from  the 
weekday  mean  KVA  line  of  Figures  7a  and  7b. 


Example  1: 

A0R-1 

CGN-38 

(AOR-1  +  CGN-38) 

^Observations 

13300  KVA 

40720  KVA 

54020  KVA 

Obser visions 

24 

24 

24  pairs 

Mean 

554.7  KVA 

1696.67  KVA 

2250.82  KVA 

554.16  +  1696.67  -  2250.82  -  2250.82  KVA 


4.3  Variance 


A  second  population  parameter  that  Is  of  Interest  in  this  study  is  i r  2  or 
the  population  variance.  Variance  is  the  measure  of  the  dispersion  of  a  popu¬ 
lation.  The  variance  within  a  population  is  the  sum  of  the  squares  of  the 
difference  between  each  particular  observation  and  the  mean.  This  sum  is  then 
divided  by  the  total  number  of  observations. 

,2  .  e  -jit  (3) 

i=l  N 

When  just  a  sample  of  information  about  the  population  is  known,  then  the 
sample  variance,  s2,  is  defined  as: 

s2  -  £  <*  ~^>2 

i-1  n  -  1  (4) 

The  n-1  is  used  so  that  the  sample  variance  will  be  an  unbiased  estimator 
of  the  population.  It  is  the  range,  of  recorded  values  that  is  of  interest  in 
PSRPM  because  it  Is  the  upper  portion  of  this  range  that  defines  the  maximum 
demand  of  electrical  power  for  an  ensemble  of  ships.  Hence,  from  this  maximum 
demand,  design  specifications  for  future  ports  and  ship  mixes  can  be  developed 
and/or  evaluated. 

4.3.1  Summation  of  Variance 

Section  4.2.1  stated  that  the  mean  of  the  sum  was  equal  to  the  sum  of  the 

means.  In  the  case  of  variance,  the  sum  of  two  variances  may  not  equal  the 


variance  of  the  entire  sample.  Example  2  uses  the  same  data  as  was  used  in 
Example  1  to  show  that  the  sum  of  two  variances  may  not  equal  the  variance  of 
the  sum. 

Example  2:  Using  Equation  3. 

Variance (A0R)  +  Variance  (CGN)  j4  Variance (A0R  +  CGN) 

275.9058  +  14375.36  *  340296.4 

14651.66  j4  340296.4 

To  permit  the  computation  of  the  variance  of  the  sum  of  demands  such  as 
the  ADR  and  CGN,  the  concept  of  covariance  shall  be  introduced. 

The  relationship  between  two  variables  may,  under  certain  conditions 
(normal  correlation)  be  described  by  way  of  a  single  number,  which  is  a 
symmetrical  function  of  the  observations.  This  single  number  is  covariance  and 
is  typically  defined  in  the  equation: 

3xy  -  —  £  <*i  -*><?!  “  30  (5) 

n-1  i=i  v  ' 

where  (xj.,  yx),  (x2,  y2),  •••»  <xn,  denote  the  observed  values  of  x  and  y. 
The  sum  of  the  products  of  the  deviations  from  the  mean  may  be  expressed  as: 

£(*i  -  x)(yt  -  y)  -  £  xiyi  - -(£*i  £  yi)  (6) 

i=l  i^l  n  i=i  i=i 

The  variance  of  the  sum  of  AOR  and  CGN  observations  can  be  expressed  as 
Equation  7. 
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"(AOR  +  CGN) 


£  I(A0R1  +  CGN1)  -  (ACTT+  CCRTl* 


i-1 


n  -  1 


(7) 


In  Equation  7  the  observed  mean  values  of  AOR  and  CGN  are  the  KVA  readings 
taken  at  common  times,  (i.e.,  0100,  0200,  0300,  ...,  2400).  The  term  (AOR  + 
CON)  Is  the  summation  of  the  means  of  the  readings  and  is  not  dependent  on  the 
time  in  which  the  readings  were  taken.  A  further  expansion  of  Equation  7 
reveals : 

,  24  _ 

s  (AOR  +  CGN)  "  21  [(AORj  ~  AOR)  +  (CGNt  ~  Ofitl)  ] 


n-1 


.^/(AORi  -  A 
i»l V  n  —  1 


-  AOR)2  +  2(AORi  -  A0R)(CGN1  -  CGN)  +  (CGNj  -  C5N)2  \ 
;  n  -  1  n-1  / 


(8) 
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Recall  that  the  variance  of  the  AOR  observations  is  21  ^A0Ri  “  A0R^ 


i=l  n-1 

and  the  variance  of  the  CGN  observations  is  ^  (CGN^  -  CGN) 

i=l  n  -  1 

The  middle  term  in  Equation  8  is  in  the  form  of  the  covariance.  Refer  to 
Equation  (5). 
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£  (A0Rt  -  A0R)(CGNi  -  COR) 


i=l 


covariance 


(9) 


n-1 


I 
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The  total  variance  of  the  bun  of  the  hourly  KVA  usages  can  be  determined 
through  the  construction  of  a  covariance  table  which  is  a  convenient  way  of 
organizing  the  variance  and  covariances.  The  present  example  illustrates  this 
fact. 


Remembering: 


2  V'  (AORj  -  SffiCV 

■  «*>  £  -  ~n  -  1 - 


2  v*  (CGNj  -  c5n)‘ 

8  (CGN)  2-  - , - 

i-1  n  -  1 


pmj  m%  (AORt  -  AOR) (CGN,  -  TO 

cov(aor,  CGN)  *  Is - 1 - * - 

i3*!  n  -  i 


The  covariance  table  is  written  as : 


®2(aor) 


C0V (CGN,  AOR) 


COV(AOR,  CGN) 


s2(CGN) 


The  COV  of  an  AOR  and  a  CGN  is,  by  its  definition,  equal  to  that  of  the 
CGN  and  AOR.  Therefore,  using  Equation  (6): 


r«w 


■  vV.  .  . 


■•*■*•’*  -  *  »  *  •  *  '  »  *  > 

»•%•*  »  '  .  "  .  c  •  *  ♦  V  ‘ 
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§2(AOR  +  CGN)  *  ®2(AOR)  +  2COV(AOR,  CGN)  +  ®2(CGN)  <10> 

which  is  the  sum  of  the  terms  in  the  above  mentioned  covariance  table. 

The  calculation  of  Equation  (6)  is  repeated  for  every  possible  combination 
of  ship  classes  In  order  to  construct  the  master  covariance  table.  Ships  are 
even  compared  to  other  ahips  in  the  same  class  to  determine  Intraclass 
covariances  (as  opposed  to  Interclass  covariances). 

A. 4  Standard  Deviation 

This  statistic  aids  in  the  description  of  the  width  of  the  distribution 
curve.  On  a  normal  curve,  1  standard  deviation,  s,  on  either  side  of  the 
distribution  mean  defines  the  curve's  inflection  points.  Between  RT  +  s  and 
x  -  s  lies  almost  68Z  of  the  observed  data,  and  between  x  +  3s  and  x  -  3s  lien 
99. 9%  of  the  data. 

The  variance  of  a  sample  grouping  of  data  carries  the  units  of  the  data 
squared  (KVA2)  and  is,  by  definition,  the  square  of  the  standard  deviation. 
The  covariance  is  difficult  to  Interpret  as  a  measure  of  the  relationship 
between  two  variables  because  its  value  is  given  in  terms  of  the  observation's 
units.  However,  dividing  the  covariance  (KVA2  in  this  study)  by  the  product  of 
the  standard  deviations,  or  square  root  of  the  product  of  the  variances  (KVA2), 
a  number  which  is  dimensionless  and  independent  of  the  units  of  measure  is 
created.  This  number  is  the  correlation  coefficient. 

4.5  Correlation  Coefficient 

Correlation  is  used  to  determine  the  degree  of  association  between  two 
items.  The  statistic  r  is  called  the  sample  correlation  coefficient.  To 
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account  for  the  fact  that  electrical  usage  generally  peaks  during  particular 
morning  and  afternoon  hours,  the  PSRPM  utilizes  the  statistical  techniques 
Involved  with  correlation.  In  all  cases  -l<r<l.  If  r  -  -1,  there  Is  a  perfect 
negative  relationship  and  if  t  ■  +1,  there  is  a  perfect  positive  relationship. 
In  the  case  of  the  PSRPM,  the  correlation  coefficient  is  a  normalised  value 
indicating  the  strength  and  character  of  the  relationship  of  electrical  power 
usage  as  a  function  of  tlae  between  any  pair  of  ships.  If  the  correlation  is 
positive  the  leplication  is  that  the  maximum  and  minimum  usages  occur  at 
approximately  the  same  time;  if  the  correlation  is  negative,  the  maximum  of  the 
usage  of  one  ship  occurs  approximately  when  that  of  the  other  is  at  a  minimum. 
Thus  a  negative  correlation  has  the  effect  of  reducing  the  variability  in  the 
sum  of  the  usages.  Figure  10  shows  the  resulting  correlation  table  for  the 
ships  present  at  Sewell's  Point  during  the  survey.  The  average  weekday  and 
weekend  correlation  coefficients  are  .75  and  .22  respectively. 

The  following  shows  an  example  of  the  computation  of  a  correlation 
coefficient.  The  ships  compared  are  once  again  the  AOR-1  and  CGN-38. 

Equation  11  is  commonly  used  to  determine  the  correlation  coefficient. 

Example  3: 


sxy  -  xy  £  y 

i=l  1*1  i*l 

24  24  24 

'  nE  KAORiKCGNj)]  ~  £  UORi)  £  (CGNj.) 
and  1=1  i=l  i=l 

n  n  2 

3xx  “  n  ~  ^  H x) 

i=l  i=l 
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Then , 


24 


24 


Q  E  UOR^)  -  (^AM^)2 

v  ■  "S/-(j&)2 

24  24 

n  £  (CGN*2)  -  (  £  CG^)2 


24 

E  A0s-< 
is1 

E  CGN, 

i-1  1 


i-1 


i-1 


13300 


24 


; (AOR,)2 
1 

I 

EXCGNj)2 
i-l 

^(AOR1)(CGNi)  -  22579225 
Sxv  -  24  (22579225)  -  ( 13300)(40720) 


40720 
24- 


7377250 

69418900 


a 


yy 


24  (7377250)  -  (13300)2 
24  (69418900)  -  (4072C)2 


325400 

164000 

7935200 


325400 


4 


-  .285  -  correlation  coefficient 


'(164000)(  7935200) 


The  circled  correlation  coefficients  in  Figure  10  show  a  computer 
generated  value  of  r  -  .23. 
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Furthermore,  th«  accuracy  of  tha  pradiccion  la  baaed  on  Che  square  of  cha 
nuaber  of  ahipa  present  in  tha  modal.  For  example,  the  nuaber  of  ships  In  the 
aodal  defines  the  dimensions  of  the  covariance  table.  N  ships  produces  a  N  x  N 
table.  The  numbers  on  the  diagonal  define  the  variances  and  the  number  on  the 
off  diagonal  are  the  covariances.  We  have: 

N  Variances,  and 
-  N  Covariances 

The  variance  of  the  sum  of  demands  would  be: 

o-^  (xj  +  X2  +  •••  +  Xq)  ■  iJV  +  (N^  “  N)twV 

where  V  is  the  average  variance  and  C6V  is  the  average  covariance. 

The  magnitude  of  the  covariances  may  be  small  but  as  the  number  of  ships, 
N,  increases  so  does  the  total  variance.  The  increase  can  be  seen  to  vary  as 
the  square  of  the  number  of  ships.  As  the  variance  increases  (as  N  increases) 
the  distribution  curve  widens  and  the  area  under  the  wings  is  increased  pre¬ 
dictability  becomes  less  precise.  See  Figure  11. 


Figure  11.  Influence  of  Covariance  on  Maximum  Resource 

Specification 
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The  value  of  Che  correlation  coefficients  not  only  gives  an  idea  of  the 

tendency  of  the  ships  to  have  corresponding  peak  and  lull  periods  of  power 

utilization,  but  they  are  indicators  of  the  predictibility  of  distribution  of 

the  maximum  demand  made  by  the  ships.  In  order  to  make  the  statement  that  r  is 

a  measure  of  the  predictibility  of  total  demand  consider  the  following: 

m  covariance 

n  2 

Sjs  <x)  3  (y) 

Assume  ship  x  and  ship  y  have  a  variance  of  1  and  a  correlation 
coefficient  of  .5.  The  covariance  would  then  equal: 


rJs2(x)  s2(y)  -  .5 


and  the  sum  of  the  terms  of  the  covariance  table  would  be: 


1  .5 


.5  1 


<rZ(x  +  y)  -  3 


If  the  correlation  coefficient  were  .1  instead  of  .5  (keeping  the 

I 

individual  variances  **  1),  the  variance  of  the  sum  would  equal  2.2  which  would 
change  the  overall  shape  of  the  final  cumulative  probability  distribution 
curve.  The  overall  shape  of  the  distribution  of  the  total  demand  for  a 
resource,  and  hence  its  predictibility,  has  a  major  effect  on  the 
specifications  of  the  capacity  of  the  resource  since  specifications  which  are 
defined  by  normal  distributions  are  normally  stated  in  a  particular  number  of 
standard  deviations  about  the  mean  (x  +  ns). 
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4.6  Conservatlveneag  Factor 

It  is  known  from  experience  that  the  resource  consumption  of  a  group  of 
ships  is  less  than  the  capacities  of  all  the  Individual  ships  added  together. 
This  is  attributed  to  both  variation  in  an  individual  ship's  demand  and  coinci¬ 
dental  loading  between  ships,  leading  to  a  certain  amount  of  diversity  in  the 
loading  between  two  or  more  ships  serviced  from  the  same  transformer.  The  DM- 
25  method  accounts  for  this  effect  by  taking  the  maximum  designed  electrical 
capacity  of  the  various  ships  and  reducing  their  sum  by  multiplier  (called  a 
diversity  factor)  dependent  upon  the  number  and  class  of  ships  serviced  by  the 
transformer. 

The  PSRPM  calculates  a  diversity  factor,  not  as  a  mechanism  for  accommoda¬ 
ting  the  diversity  effect  as  with  the  use  of  the  DM-25  diversity  factor,  but 
rather  as  an  illustrative  parameter.  The  actual  capacity  itself  is  predicted 
by  the  PSRPM  and  the  diversity  factor  illustrates  how  the  capacity  for  a  ship 
group  relates  to  the  sum  of  the  capacities  for  the  individual  ships. 

The  diversity  factor  depends  strongly  on  the  criterion  established  fot 
relating  the  capacity  to  the  parameters  of  the  demand  distribution.  This  was 
accomplished  through  a  number  defined  as  the  conservativeness  factor  (K)  which 
is  the  number  of  standard  deviations  above  the  mean  demand  at  which  the  capaci¬ 
ty  requirement  is  set.  The  larger  the  K  factor,  the  more  conservative  the 
estimate  of  required  capacity. 

The  evolution  of  the  conservativeness  factors  for  the  PSRPM  may  be  found 
in  greater  detail  in  A  Comparison  of  Two  Methods  for  Estimating  the  Electrical 
Supply  Capacity  of_  Navy  Piers,  1/83.  As  a  result  of  this  study,  the  recommen¬ 
dation  was  made  to  continue  the  DM-25  process  as  a  direct,  simple,  and  conser¬ 


vative  process  for  designing  the  power  capacity  of  permanent  installations. 
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V.  PSRPM  -  PHASE  II  COMPUTER  PROGRAMS 

5.0  Introduction 

Phase  I  of  the  PSRPM  utilized  a  suite  of  five  data  base  programs  and  two 
data  analysis  programs  in  addition  to  computer  aided  graphics  package  to 
generate  the  probability  (Demand  >  Capacity)  versus  Capacity  curves.  Phase  II 
reduced  the  number  of  data  base  programs  to  one,  and  modified  one  of  the  data 
analysis  programs  to  permit  future  refinements. 
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5.1  PSP2  FORTRAN  Programs 

The  program  PSP2  is  an  enhanced  version  of  program  PSP1.  A  listing  of 
PSP1  may  be  found  in  the  Summary  Report;  An  Application  of  the  Port  Systems 
Requirements  Prediction  Methodology  (PSRPM)  Using  Electrical  Measurements 
Recorded  at  Sewell's  Point  (12-76  to  1-77).  PSP2  generates  realizations 
similar  to  PSP1  but  substitutes  electrical  usage  data  obtained  from  the 
observations  made  at  Sewell'3  Point  in  place  of  DM-25's  Present  Demand  values 
in  calculating  the  total  KVAs  used.  There  is  one  data  base  program  and  one 
data  analysis  program  which  assist  PSP2.  Complete  listings  of  the  programs  and 
sample,  of  the  data  files  are  located  in  Appendix  A. 

5.1.1  Data  Base  Program 

The  surveyed  data  was  formatted  for  FORTRAN  and  entered  into  a  program 
called  KVADB.  This  program  created  the  data  base  KVA.DBS  (KVA  Data  Base)  which 
contains  the  ship  class,  number  of  days  in  port,  date,  daily  hourly  KVA 
readings,  and  mean  hourly  KVA  usage,  and  summed  mean  hourly  KVA  usage.  The 
data  base  distinguishes  weekday  from  weekends  by  the  use  of  a  minus  sign  in 
column  1.  Appendix  A  contains  a  printout  of  KVADB  and  KVA.DBS,  and  Figure  12 
presents  a  schematic  of  this  logic. 
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Figure  12.  Data  Base  Logic  Flow 
5.1.2  Data  Analysis  Program 

A  program  separate  from  PSP2,  named  PSP0  was  created  to  construct  the 
covariance  and  correlation  tables.  PSP0  was  made  separate  from  PSP2  to  aid  in 
the  incoropratlon  of  refinements  to  the  statistical  techniques  which  were  used 
to  calculate  values  in  each  table.  Refinements  and  changes  to  procedures  can 
simply  be  made  to  PSP0  without  having  to  modify  the  master  PSP2  program. 

Figure  13  provides  a  view  of  the  computer  process.  PSP0  uses  data  from 
KVA.DBS  to  generate  the  Master  Covariance  and  Correlation  Tables.  Information 
from  the  KVA.DBS  and  these  tables  is  inputted  to  PSP2  which  generates  the  data 
points,  from  the  realizations,  for  the  P(D>C)  vs.  C  curves  for  different  port 
scenarios. 
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6.0  Introduction 

All  of  the  statistics  discussed  in  Section  IV  have  a  relevance  to  the 
prediction  of  the  overall  demand.  The  mean  for  each  class  of  ships  defines  the 
center  of  the  distribution  for  that  class  and  the  total  variance  defines  the 
shape  of  the  distribution  in  the  form  of  the  standard  deviation.  The 
correlation  coefficient  provides  an  Indication  of  the  degree  of  association 
between  the  peak  electrical  demand  periods  of  any  pair  of  ships. 

6.1  Objective 

The  objective  of  this  section  is  to  discuss  the  application  of  the 
statistics  presented  in  Section  IV  and  to  present  the  importance  of  the 
contribution  of  the  covariance  in  the  shaping  of  the  distribution  curve. 


6.2  Covariance 

Hourly  KVA  readings  were  taken  over  the  period  of  several  days  for  28 
ships  in  22  classes  at  Sewell's  Point.  This  data  was  used  in  this  study  to 
demonstrate  a  methodology  to  predict  future  port  utilization  requirements. 
PSRPM  is  a  demonstration  of  a  particular  type  of  prediction  methodology  which 
utilizes  the  means  and  variances  associated  with  the  hour  to  hour  and  day  to 
day  electric  power  consumption  utilizations.  The  computation  of  the  sum  of  the 
variances  which  define  the  future  requirements  is  not  a  trivial  task.  In  order 
to  simplify  the  computation  a  physical  manipulation  of  the  terms  in  the 
quadradlc  equation  (see  Equation  7}  permits  the  summation  of  the  variances 
and  covariances  belonging  to  two  sets  of  ship  data.  Executing  this  procedure 
for  every  pair  of  ships  (784  Phase  II  combinations)  results  in  a  28  x  28 
covariance  table.  The  Master  Covariance  Table  is  presented  in  Figure  14. 
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Th«  Master  Covariance  Table  is  unique  because  it  provides  a  visual  and 
easily  accessible  representation  of  the  variation  in  KVA  utilisation  of  any 
particular  class  nix  of  ships. 

Take  for  example  a  group  of  ships  made  up  of  S2,  S^,  S5,  and  Sm.  The 
value  of  the  ensemble  variance  would  equal  the  sum  of  the  corresponding 
variance  and  covariance  terms  taken  from  the  Master  Covariance  Table.  Figure 
IS  exhibits  this  example. 
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Figure  15.  Use  of  Master  Covariance  Table 


Eusemble  Variance  would  then  equal  V_  „  +  V,  ,  +  V,.  +  V  + 

2,2  4,4  5,5  m,m 

2  COV  ,  +  2  COV  ,  +  2  COV  +2  COV.  ,  +  2  COV.  +  2  COV.  which  is 
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6.3  Correlation  Coefficient 


If  th«  correlation  between  the  pairs  of  ships  analysed  were  sero  then  by 
definition  (see  Equation  12)  there  would  not  be  a  covariance  term  in  the 
quadradic  equation  defining  the  sum  of  the  variances  (see  Equation  9).  The 
absence  of  covariance  would  permit  the  summation  of  the  individual  variances  in 
defining  the  shape  of  the  distribution  curve.  Section  IV  showed  that  there  is 
in  fact  correlation  between  ship's  time  of  resource  usage  and  thus  the 
contribution  of  the  covariance  term  to  the  shape  of  the  distribution  is 
important.  Figure  11  presents  a  representation  of  how  the  maximum  demand  can 
be  Influenced  by  the  incorporation  of  covariance. 


6.4  Determination  of  Variations 

The  major  portion  of  the  PSRPM  work  was  in  the  identification  and 
determination  of  the  variances  in  K VA,  utilization  that  occurred  between  ships 
at  cold  iron.  The  data  sheets  which  were  constructed  for  the  survey,  see 
Figure  7a  and  7b,  permitted  the  calculation  of  the  hourly  and  daily  means  and 
variances.  Figure  16  is  presented  as  an  aid  in  che  discussion  of  the 
determination  of  the  total  variance  in  a  ship's  power  utilization. 
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Figure  16.  Notional  Sewell's  Point  Data  Sheet 


From  the  hourly  means,  xh,  a  mean,  xhd  and  a  variance  s2h>  can  be 
calculated  which  represents  the  distribution  of  the  power  utilization  over  the 
period  of  one  day.  The  daily  means,  x^[,  produce  a  mean,  x^  which  is  the  same 
as  the  mean  of  'he  hourly  means, "x^i  (see  Section  4.2.1).  But  the  variance 
that  is  associated  with  the  daily  means,  s2d,  is  not  equivalent  tu  S2h»  It  is 
at  this  point  where  an  lndepth  look  at  the  variances  is  warranted. 

A  simplistic  approach  to  the  development  of  the  methodology  would  be  to 
neglect  any  correlation  between  the  ships'  power  utilization  and  calculate  the 
mean  and  sum  the  standard  deviations  for  each  ship  class.  These  two  statistics 
would  define  the  KVA  utilization  distributions  for  the  various  classes  of 
ships.  An  assumption  of  a  zero  correlation  would  result  in  the  values  of  the 
off  diagonal  elements  of  the  covariance  table  being  zero.  Thus  the  total 
variance  of  the  distribution  would  result  from  the  sum  of  the  diagonal 
elements.  The  value  of  the  standard  deviation  in  a  N  x  N  covariance  table 
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would  be  influenced  by  the  factor  N.  Thi*  procedure  however  would  deny  to  the 
study  tha  lnfluanca  of  tha  operational  and  environmental  variations. 

Tha  nose  consarvativa  approach,  which  would  maximise  tha  lnfluanca  of  tha 
operational  and  environmental  variations,  aakas  tha  assuaption  that  thara  is  a 
parfact  corralation  (r»l)  batwaan  tha  powar  utilisations  of  tha  various  ship 
classas.  This  assuaption  craatas  a  valua  of  1  for  tha  diagonal  and  off 
diagonal  aleaants  of  tha  covarianca  tabla.  Hence,  tha  valua  cf  tha  total 
variance  would  propogata  as  the  square  of  tha  nuabar  of  ship  classas  analysed. 
(N2)  *nd  tha  standard  deviation  would  propagate  as  tha  sum  of  tha  standard 
deviations.  Tha  resulting  statistics  which  dafina  tha  utilisation  distribu¬ 
tions  would  consist  of  a  class  naan  and  a  standard  deviation  influenced  by  +  N. 
Exanpla  4  Illustrates  the  bound ry  conditions  established  by  0  <  t  <  1. 

Example  4: 

N  -  10  classes  of  ships  each  with: 

s  -  100  KVA 

s2  -  100  KVA 

s  -  10  KVA 

With  no  corralation  (r  -  0)  the  summation  of  the  10  diagonal  elements  of 
the  covariance  table  results  in  a  total  variance  of: 

10  (100)  -  1000  -  s2 

thus 

s  ■  31.6 
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With  perfect  ensemble  correlation  the  summation  of  all  the  elements  of  the 
covariance  table  results  in  a  total  variance  of: 

s2  -  (100)  (100) 

thus 

s  -  100 

Both  conditions  would  have  the  sane  naan  but  different  standard  devia¬ 
tions. 

6.4.1  Variations  in  KVA  Utilisation 

Two  factors  nay  contribute  to  the  variations  in  ship's  electrical  power 
consumption.  Thane  factors  are  the  operational  conditions  of  the  ship  and  the 
envlronaental  conditions  present  at  the  port.  Variations  in  KVA  utilisation 
occurred  froa  weekday  to  weekend.  The  operations  profiles  for  the  majority  of 
the  ships  were  based  on  a  Monday  to  Friday  work  week.  Noticeable  reductions  in 
the  power  requirements  of  each  ship  occurred  on  weekends  primarily  due  to  the 
crew  not  conducting  routine  training  and  maintenance  operations  on  Saturdays 
and  Sundays.  The  differences  between  weekday  and  weekend  power  requirements 
were  great  enough  to  require  a  segregation  of  the  data.  Only  weekday  data  was 
used  in  the  construction  of  covariance  and  correlation  tables. 

6.4.2  Operational  Variations 

Figure  17  shows  the  probability  of  incurring  the  maximum  daily  demand  as  a 
function  of  the  time  of  day.  It  Illustrates  that  peak  usages  usually  occur  at 
1000  hours  and  at  1400  hours.  This  is  most  likely  due  to  the  preparation  of 
meals,  training,  and  the  required  system  check  outs.  The  correlation  that 
occurs  between  classes  of  ships  froa  hour-to-hour  is  incidental.  It  is  an 
explicit  relationship,  strongly  dictated  by  present  managerial  practices. 
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6.4.3  Environmental  Variations 

Environmental  factors  such  as  temperature,  wind  speed,  and  cloud  cover 
effect  the  power  requirements  of  a  mix  of  ships  on  a  day  to  day  basis.  Changes 
in  these  conditions  do  not  usually  happen  rapidly  enough  to  be  classified  as 
time  of  day  variations  so  they  are  considered  to  be  day  to  day  variations.  The 
influence  of  the  environment  can  be  pictured  as  an  addition  to  or  subtraction 
from  the  ship's  required  operational  loads.  That  is,  the  ships  will  carry  out 
their  normal  day  to  day  activities  and  will  turn  the  heat  up  on  cold  days  or 
use  the  air  conditioner  on  warm  days,  thus  increasing  the  total  overall  power 
utilization.  The  effect  of  the  environment  is  the  largest  cc  ?.r.ributing  factor 
to  the  total  variance  of  the  distribution. 

6.5  Problems  in  the  Data  Sets 

Hourly  readings  of  electrical  power  utilization  were  recorded  for  the 
ships  present  at  Sewell's  Point  from  20  December,  1976  to  21  January,  1977. 
Numerous  classes  of  ships  were  present  in  port,  but  not  all  at  the  same  time, 
and  their  time  spent  at  cold  iron  was  not  necessarily  continuous  from  day  to 
day.  Problem  areas  in  the  data  sets  occurred  in  four  areas:  (1)  coinciden¬ 
tally  of  data,  (2)  discontinuity  of  data,  (3)  weekday  to  weekend  variations, 
and  (4)  cold  iron  status.  The  only  method  available  to  cope  with  inconsistent 
or  discontinuous  data  was  for  the  investigator  to  make  educated  guesses  on  what 
the  missing  data  would  have  been  and  on  what  effects  missing  data  had  on  the 
final  results.  An  elaboration  on  each  of  the  problem  areas  is  included  below. 

6.5.1  Coincidence  of  Data 

Each  ship  maintained  its  individual  schedule  of  operations  and  arrived  and 
departed  according  to  operations  orders.  Figure  18  shows  a  timeline  of  the 
presence  of  each  ship  for  which  data  was  recorded.  The  lack  of  overlapping 
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data  prevented  the  complete  incorporation  of  the  day-to-day  or  environmental 
variations  into  the  covariance  table. 


6.5.2  Discontinuity  of  Data 

Interruptions  in  data  sets  occurred  periodically.  The  two  major  reasons 


for  discontinuous  data  sets  were  the  failure  of  data  recorders  and  the  shifting 
of  the  ship '8  desire  to  use  ship's  power  versus  port  power. 


6.5.3  Cold  Iron  Status 

The  entire  PSRPM  study  has  assumed  that  if  a  ship  was  receiving  power  from 
the  port  than  it  was  at  cold  iron.  Due  to  the  location  of  the  electrical 
measuring  equipment  (at  the  transformers),  it  was  not  possible  to  determine  if 
the  ships  drawing  port  supplied  power  were  in  fact  at  cold  iron.  The  knowledge 
of  whether  a  ship  is  at  cold  iron  or  not  has  an  impact  on  the  calculations  and 
assumptions  leading  to  the  determination  of  the  maximum  demand.  If  ships  are 
using  supplied  power  and  their  probability  of  usage  has  not  been  taken  into 
consideration,  then  the  estimates  arrived  herein  are  low. 


6.6  Summary 


The  statistics  discussed  in  this  report  each  had  an  impact  on  the  total 


electrical  demand  of  the  ensemble.  The  mean  established  the  central  region  of 
utilization,  the  correlation  coefficient  showed  that  there  was  more  than 


individual  variance  contributing  to  the  distribution  of  curve  defining  the 


power  utilization,  and  the  covariance  was  the  contributing  factor  to  the 


increase  in  total  variance 
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VIZ.  CONCLUSIONS 

7.0  Introduction 

The  primary  output  of  Phase  II  is  similar  in  format  to  that  of  Phase  I.  A 
graph  displaying  the  probability  that  the  demand  for  electric  power  will  exceed 
capacity  versus  capacity  was  generated.  Three  scenarios  are  displayed  in 
Figure  19.  Scenario  1  is  basically  a  Phase  I/PSP1  output  using  the  ships  that 
were  present  at  Sewell's  Point  on  22  December  1976.  Scenario  2  used  PSP2  and 
only  the  mean  KVA  utilizations  as  Inputs  and  Scenario  3  included  is  representa¬ 
tive  of  the  curve  if  the  Master  Covariance  Table  values  were  incorporated  into 
PSP2. 


7.1  Objective 

The  object  of  this  section  is  to  present  the  conclusions  reached  as  a 
result  of  the  PSRPM  Phase  II. 

7.2  Scenario  One  Sewell 'a  Point  -  Using  DM- 2 5  Values 

A  principle  output  of  the  PSRPM  is  an  estimate  of  the  probability  that  a 
group  of  ships  will  exceed  the  electrical  supply  capability  of  the  port 
facility.  This  concept  was  demonstrated  in  Phase  I  using  Present  Demand 
specifications  stated  in  DM-2S.  In  Scenario  One  PSP1  was  run  using  the  ships 
that  were  present  at  Sewell's  Point  on  22  December  1976  as  a  data  base  and 
their  respective  DM-25  values  as  KVA  utilizations.  The  distributions  created 
by  the  realizations  involving  the  22  December  ships  is  represented  in  Figure  19 
as  dots. 
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7.3  Scenario  Two  Sewell's  Point  -  Using  Mean  KVA  Values 

A  second  sec  of  values  were  calculaCed  and  plotted  which  reflect  the 
ships'  utilization  of  the  calculated  mean  KVAs.  The  mean  values  for  each  ship 
class  were  calculated  over  the  surveyed  tine  period*  Phase  II  generated  a 
series  of  random  numbers  to  determine  the  number  of  ships  present  in  each  class 
per  realization  in  accordance  with  the  Binomial  Theorem.  Next,  the  mean  KVA 
values  corresponding  to  each  ship  class  were  multiplied  by  the  number  of  ships 
in  that  class  present  in  the  realization  and  summed.  Like  PSP1,  PSP2  executed 
100  realizations  and  listed  the  total  KVA  values  in  decending  order.  The 
results  of  their  distribution  and  the  values  corresponding  to  the  Central  Limit 
Theorem  are  plotted  in  Figure  19  as  circles  (o)  and  a  solid  line  respectively. 

7.4  Scenario  Three  Sewell's  Point  Using  Covariance  Table  Values 

The  third  curve  on  the  graph  represents  what  the  probability  of  failure  to 
meet  demand  curve  would  look  like  if  all  of  the  factors  contributing  to  varia¬ 
tion  were  incorporated  into  the  covariance  values.  Since  the  variations  due  to 
the  environment  were  not  Included  the  curve  was  shifted  slightly  to  the  right 
from  where  it  would  actually  be  located.  This  compensates  for  not  incorpora¬ 
ting  this  variance  in  the  model.  PSP2  generated  a  series  of  random  numbers 
to  determine  how  many  ships  of  each  class  are  to  be  present  in  each  realiza¬ 
tion.  To  determine  the  KVA  demand  of  the  ships  present,  PSP2  is  formatted  to 
initiate  a  slightly  more  advanced  procedure  based  on  the  distributions  created 
by  each  ship's  mean  and  the  single  total  variance  which  was  derived  through  the 
summation  of  the  elements  in  the  covariance  table  which  are  particular  to  that 
realization.  If,  in  the  process  of  constructing  a  realization,  there  was  not  a 
Sewell's  Point  calculated  covariance  available  PSP2  substituted  representative 
statistical  values  for  the  covariance.  These  representative  values  were  deter- 
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mined  by  using  a  mean  equal  to  60%  of  the  DM-25  value  and  a  standard  deviation 
of  6%  of  the  DM-25  value  in  the  computation  of  the  covariance.  A  correlation 
coefficient  of  .8  was  used  as  the  final  term  in  the  calculation  of  the  repre¬ 
sentative  covariance.  These  particular  values  for  x,  s,  and  r  were  determined 
through  analysis  and  comparison  of  the  Sewell's  Point  data  to  DM-25  values. 
These  points  on  the  cumulative  distribution  curve  are  plotted  as  (*)  in  Figure 
19. 


7.5  Interpretation  of  Scenario  Curves 

Each  scenario  curve  can  be  analyzed  separately  to  reveal  particular 
characteristics  of  demand  specification  versus  probability.  Or  all  three  of 
the  curves  can  be  analyzed  as  a  whole  to  present  a  picture  of  how  the  PSRPM 
predicted  values  compare  to  the  DM-25  predicted  values,  and  how  the  addition  of 
covariance  changes  the  shape  of  the  curves,  particularly  in  the  distribution 
wing. 

7.5.1  PSP2  Mean  Curve  Compared  to  PSP2  DM-25  Curve 

It  has  always  been  surmised  that  DM-25  over  estimates  a  port's  electrical 
power  requirements.  But  the  safety  margin  built  into  DM-25  associated  with 
this  over  designing  has  never  been  determined.  Figure  20  presents  KVA  values 
of  the  PSP2  Mean  Curve  and  the  PSP2  DM-25  curve  for  constant  probabilities  of 
failure  to  meet  demand.  From  these  values  it  appears  as  though  the  DM-25 
process  incorporates  a  safety  factor  of  approximately  2. 
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7.5.2  PSP2  Mean  Curve  Compared  to  PSP2  Covariance  Included  Curve 

The  PSP2  Mean  Curve  compared  Co  Che  PSP2  Covariance  Included  curve  shows 
Che  influence  of  Che  operacional  variances  only.  As  stated  before,  the 
environmental  variations  were  not  incorporable  because  of  non-coincident  data. 
The  slight  difference  in  the  mean  hourly  KVA  utilisations  from  weekdays  to 


Acceptable  Failure  to 

PSRPM 

DM- 25 

Safety 

Meet  Demand  Rate 

KVA 

Requirement 

KVA 

Requirement 

Factor 

0 

9700 

21700 

2.23 

.1 

8500 

17500 

2.06 

.2 

7200 

15000 

2.08 

.4 

6000 

12600 

2.10 

.6 

4700 

11500 

2.45 

.8 

3600 

7500 

2.08 

Figure  20.  Determination  of  Safety  Factor 


weekends  (7.3%  greater  on  weekdays  for  AOR  (38  KVA),  and  1.25Z  greater  on 
weekdays  for  CGN  (21  KVA)  reveal  that  as  majority  of  the  power  utilized  is  in 
fact  used  for  environmental  and  hotel  purposes.  As  stated  in  Section  7.4  the 
Covariance  Included  Curve  was  shifted  slightly  to  the  right  to  compensate  the 
study  for  not  Incorporating  the  dally  environmental  variations.  This  shift  was 
justified  since  the  output  of  the  Master  Correlation  Table  shows  that  there  is 
a  very  high  correlation  between  most  of  the  classes  of  ships  even  without 
considering  the  environmental  variations.  Keeping  this  in  mind,  the  Scenario  3 
curve  tends  to  approach  the  maximum  staudard  deviation  boundary  condition,  r-i. 
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Studies  and  surveys  to  date  have  been  mainly  concerned  with  the  mean  power 
consumed  by  each  single  ship  in  port.  Periodically  maximum  demands  were 
recorded  but  no  correlation  to  the  time  of  the  occurrence  of  the  maximum 
demands  has  been  made.  The  PSRFM  has  tried  to  "statistically  blend"  the  means 
and  variances,  into  a  distribution  which  is  different  from  a  distribution  based 
solely  on  the  mean  utilisations.  The  resulting  widening  of  the  wing  redefines 
the  maximum  power  specifications  of  a  port  given  a  particular  mix  of  ships. 

7.6  PSRPM  Versus  DM-25 

Even  though  the  PSRPM  process  provides  an  accurate  prediction  capability  » 
it  is  a  complex  procedure  to  use.  Additionally,  It  makes  no  allowance  for  the 
increase  in  demand  that  occurs  in  the  life  span  of  a  typical  Navy  ship.  PSRPM 
is  useful  and  economical  when  an  estimate  of  actual  utilization  would  meet  the 
needs  of  the  circumstances  being  investigated*  The  methodology  can  aggregate 
individual  ship  power  usage  data  to  predict  the  usage  by  various  combinations 
of  ships  and  thereby  is  capable  of  predicting  power  use  at  the  pier  level. 

The  comparison  of  a  PSRPM  prediction  compared  to  a  DM-25  prediction 
results  in  an  overall  reduction  in  the  predicted  power  demand.  Thus,  a  design 
capacity  derived  from  DM-25  appears  to  be  conservtive. 


r. 


7.7  Limitations  to  the  PSRPM 

The  PSRPM  output  provides  a  demonstration  as  to  the  potential  of  the 
statistical  techniques  available  for  port  system  analysis.  The  statistical 
techniques  used  in  the  development  of  this  methodology  incorporate  the 
correlation  concept  which  in  turn  allow  for  the  computation  of  statistically 
based  diversity  factors  and  an  estimation  of  the  probability  of  failure  to  meet 
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demand  given  the  covariance  table  values.  Although  the  data  was  sufficient  to 
refine  PSP1  into  PSP2  the  methodology  is  limited  in  several  areas. 

The  Sewell's  Point  Survey  measured  power  utilization  at  the  transformer 
level.  Without  the  supporting  details  of  usage  conditions  (cold  iron  status, 
stand  down  time,  testing  of  ship's  systems,  general  overhaul,  temperatures  and 
lack  of  sufficient  overlapping  data,  etc.)  the  data  reduction  is  subject  and 
limited  to  broad  statistical  interpretation.  Since  the  lata  was  collected  over 
the  Christmas  holiday  season,  when  many  personnel  were  on  leave  and  ship 
activity  was  low,  the  power  consumption  cannot  be  considered  representative  of 
the  yearly  normal.  In  addition,  it  is  uncertain  whether  an  hourly  sample 
period  provides  sufficient  detail  to  identify  peak  loads  of  short  duration. 

The  applicability  cf  the  Binominal  and  Central  Limit  Theorems  impose 
limitations  to  the  study.  An  understanding  of  the  operational  profiles  of  the 
ships,  such  as  the  time  in  port  as  compared  to  deployment  time  is  essential  in 
the  assumption  of  a  binominal  distribution.  Further  investigation  on  the 
Influence  of  one  ship  s  electrical  consumption  affecting  another  ship's 
consumption  and  on  the  accuracy  of  the  assumption  that  if  a  ship  is  not  its 
houeport  then  its  presence  is  substituted  with  a  similar  ship  at  the  homeport, 
needs  to  be  accomplished.  The  Central  Limit  Theorem  made  the  assumption  that 
the  data  could  be  represented  by  a  normal  distribution.  The  sample 
distributions  obtained  in  this  study  are  generally  closer  to  the  normal 
distribution  in  the  median  region  with  inaccuracies  in  the  wings. 

Taking  into  consideration  the  limitations  stated,  and  realizing  that  the 
PSRPM  is  the  demonstration  of  a  concept  which  is  intuitively  correct  with 
problems  to  be  resolved,  it  is  concluded  that  the  results  of  this  methodology 
be  considered  preliminary  and  that  their  use  for  planning  estimations  be  sub¬ 
ject  to  this  fact. 
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VIII.  RECOMMENDATIONS 


8.0  Discussion 

Th«  Pore  Systems  Requirements  Prediction  Methodology  waa  developed  to 
demonstrate  the  ability  to  make  predictiona  about  port  utility  aervlce 
requirements.  The  application  of  the  techniques  developed  in  the  PSRPM  is  not 
limited  to  the  electrical  power  resource  but  la  applicable  to  every  aervlce 
reeource  the  port  has  to  deliver  separately  or  in  combination.  A  foundation 
for  the  future  development  of  prediction  methodologies  was  made  through  the 
lessons  learned  in  the  overall  development  of  the  PSRPM.  Based  on  PSRPM. 
enhancements  to  the  decision  making  process  and  to  the  design  of  future  data 
gathering  programs  are  envisioned  in  experimental  design;  understanding  of 
present  design  methods;  and  economic  analysis. 

Phase  I  was  the  initial  attempt  in  developing  a  prediction  methodology. 
It  validated  the  theoretical  basis  by  comparing  calculated  data  to  the  Central 
Limit  Theorem  Data.  The  calculated  data  was  very  clean  in  nature  since  it  was 
composed  entirely  of  theoretical  means.  The  output  of  Phase  I  proved  to  be 
usable  and  quite  understandable  in  its  application  to  the  objective  of  the 
study.  Phase  II  built  on  the  first  phase  by  incoropratlng  actual  measured  data 
and  Introducing  hourly  and  daily  utilization  variances.  The  calculated  values 
once  again  tracked  very  closely  with  the  theoretical  Central  Limit  Theorem 
values.  An  understanding  of  ship  correlation,  covariance,  and  the  influence  of 
large  numbers  of  ships  on  the  overall  outcomes  of  a  prediction  has  now  been 
developed.  Phase  II  presents  areas  where  Improvements  to  both  the  experimental 
process  and  the  application  of  statistics  may  be  improved.  The  next  logical 
step  is  the  development  of  a  more  refined,  accurate  prediction  methodology, 
possibly  a  Phase  III,  which  would  involve  a  validation  by  direct  comparison  of 
its  conclusions  with  monitored  parameters. 
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8.1  Objective 

The  objective  of  this  section  is  to  present  recommended  applications  of  a 
PSRPM  basad  procadura. 
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8.2  Bxporiswntal  Design 

Tha  PSRPM  product  was  davalopad  to  incorporata  tha  Sewell's  Point  data. 
Tha  statistical  tachniques  usad  wara  adapted  to  what  information  was  available, 
and  tha  results  are  a  product  of  this  adaptation.  Tha  specific  area  whan  tha 
PSRPM  has  contributed  to  future  experimental  designs  is  in  tha  planning  of  data 
acquisition  process.  Sub  areas  of  tha  data  acquisition  process,  such  as  whan 
to  take  tha  data,  what  data  to  taka,  and  how  many  ships  should  be  instrumented, 
have  all  been  addressed  in  the  development  of  the  PSRPM. 

8.2.1  When  to  Record  Data 

The  Sewell's  Point  data  showed  that  generally,  the  maximum  power 
requirements  occurred  in  the  late  morning  and  in  the  early  to  mid  afternoon. 
Knowing  this,  it  does  not  appear  necessary  to  record  data  around  the  clock. 
Instead,  measurements  of  more  frequent  Intervals  or  of  a  continuous  nature 
taken  during  the  observed  peak  usage  periods,  would  display  short  periods  of 
maximum  power  utilization.  Since  the  Sewell's  Point  data  was  acquired  in  the 
wintertime,  the  complete  spectrum  of  the  ships'  power  requirements  for  heating 
and  cooling  has  not  been  completed.  It  appears  feasible  to  schedule  the  next 
data  gathering  exercise  during  the  season  which  traditionally  has  the  greatest 
Influence  on  the  base's  overall  power  requirements.  Sampling  rate  requirements 
have  been  reviewed  and  it  has  been  determined  that  operationally  based 
variations  require  more  rapid  sampling  than  the  hourly  case.  Environmentally 
Induced  variations  can  generally  be  estimated  and  should  be  sampled  less 
frequently. 
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8,2*2  What  Data  to  Record 

An  in  dapth  analysis  of  tha  correlation  coafficisnta  and  tha  naan  power 
utilisations  of  tha  various  classes  of  ships  will  aid  the  experiment  designer 
in  developing  a  more  efficient  scheae  of  ship  instrumentation.  However,  ship 
level  power  utilisation  information  is  required.  The  next  logical  step  would 
be  to  measure  the  utilisation  at  tha  ship  level,  perhaps  at  the  ship's  buss, 
to  determine  if  the  ship  was  in  fact  at  cold  iron.  This  would  permit  the 
extrapolation  of  the  data  to  all  of  the  ships  of  a  similar  type,  and  hence, 
port  scenarios  can  be  synthesised  and  analysed, 

8,3  Assistance  to  Port  Master  Planners 

The  concept  of  forecasting  the  future  can  be  seen  to  be  very  complex.  One 
can  only  make  educated  guesses  on  technological  advancements  end  requirements 
of  the  future.  The  PSRPM  provides  understandable  and  almost  immediate  feedback 
on  the  "what  if's"  of  the  future.  Additionally,  with  refinement,  an 
interactive  capability  is  envisioned.  This  study  can  benefit  the  Port  Systems 
Project  by  providing  analytical  procedures  for  utility  demand  data  acquisition 
projects.  With  futher  refinement,  PSRPM  could  be  used  to  assess  the  impact  of 
selected  MILCON  Projects  and/or  serve  as  a  basis  of  an  acquisition  strategy. 
Finally,  as  familiarity  with  the  process  is  gained,  procedures  could  be 
developed  which  would  present  NAVFAC  and  fleet  planners  with  an  analytical 
review  of  existing  and  planned  facilities  on  an  interactive  basis. 

t  -‘■N 

8,3.1  Present  Design  Procedures  -  DM-25 

The  present  process  for  establishing  the  electrical  power  service 
requirements  for  a  port  involves  the  blending  of  knowledge  about  a  berthing 
plan  with  the  DM-25  process.  Until  now,  it  was  not  clear  to  port  planners  how 
conservative  DM-25  was.  Results  from  the  PSRPM  show  that  DM-25  designs  with  a 
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safety  margin  of  approximately  2.  The  P3RPM  can  provide  some  insight  to  how 
large  the  power  distribution  system  should  be  without  approaching  an 
overdesigned  specification.  Because  of  its  simplicity,  directness,  and  built 
in  conservativeness  the  DM-25  procedure  should  remain  the  primary  design 
approach. 


8.3.2  Economic  Analysis 

Capital  investment  strategies  for  port  electrical  utilities,  time  value  of 
money,  and  rate  of  return  are  all  based  on  two  questions;  how  much  power  is 
required  in  the  future  and  for  how  long  is  this  present  supply  level  adequate? 
The  PSRPM  can  provide  assistance  in  estimating  how  long  a  particular  power 
system  can  deliver  the  maximum  load  specification.  Given  an  acceptable  failure 
to  meet  demand  rate  and  the  known  mix  of  ships  a  time  of  useful  life  of  the 
present  power  system  can  be  established.  From  this  information,  the  most 
economical  use  of  funds  available  can  be  determined. 

8.3.3  Lead  Time  Estimation 

When  the  decision  is  made  to  modify  the  port's  existing  facilities,  the 
PSRPM  can  provide  an  estimate  to  the  planners  as  to  how  long  the  present  system 
can  satisfactorily  supply  power.  This  information  gives  the  port  planner 
guidance  on  when  to  schedule  the  required  construction  projects. 

8.4  Summary 

The  PSRPM  project  demonstrated  that  the  statistical  techniques  developed 
can  allow  port  planners  to  establish,  with  some  confidence,  the  utility  usage 
demanded  for  various  combinations  of  ships  which  are  expected  to  be  served.  A 
strong  relationship  between  the  data  acquisition  process  and  the  validity  of 
the  PSRPM  can  be  implied. 
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PHASE  II  COMPUTER  PROGRAMS 
LISTINGS  AND  EXAMPLES 
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K7AD3 

KVADB  is  a  FORTRAN  program  which  creates  the  data  base  KVA. DBS.  It 
contains  the  ship  class,  number  of  days  in  port,  date,  mean  hourly  KVA  usage, 
and  summed  mean  hourly  KVA  usage. 


.DIR  KVADB ♦ FOR 


i 

i 

m  » 

i  ;0 

■  .1 


XDRTLKE  Non -ex  i  stont  file  KVADB .  FOR 
yDIR  CKV\VN\5734 » 1303  KVADB . FOR 

KVADB  FOR  4  <057>  29-Dec-81  NOAFt  C5734»1301 

♦  TV  C5734  *  1303  KVADB .FOR 
PROGRAM  KVADB 
IMPLICIT  INTEGER  (A-2) 

REAL  FACTOR 

DIMENSION  NAME2 ( 2 > ? D < 30 y 25 ) t M ( 2 r 25 ) » S ( 2 » 25 ) 

.10  TYPE  100 

READ ( 5 t 1 10  r ERR- 1 0 )  NAME 

OPEN ( UNIT-6 1 ACCESS” ' 3EQIN ' y  FILE-NAME , DEVICE- ' DSN ' ) 

OPEN ( UNIT-7 v  ACCESS- ' APPEND ' ,FI LE~ ' KVA . DBS ' r DEVICE- ' DSK ' ) 
READ (6k 120)  NAMES  *  DAYS 
DO  20  I  :k1»  2 

READ ( 6 f 1 30)  < M < I ?  J ) y J“1 . 25 > 

20  CONTINUE 

DO  30  1*1 » DAYS 

READ ( 6 1 130 )  ( D ( I » J ) » J= 1 y  25  > 

30  CONI INUE 

C  DETERMINE  CONVERSION  FACTOR 
iil*0 
IM1-0 

DO  40  J  ~  1  ?  D  A  Y  8 

IF  (D(jyl)  ♦ EQ ,  -1>  GOTO  40 
Ml  “Ml  +D  <  J  .*  2  ) 

IM1-IM1 +1 
40  CONTINUE 

Ml-IFIX  ( FLOAT  ( M:l )  /F  LOAT <  IM1 )  ) 

FACTO REFLOAT ( M (1 y  2 ) ) /FLOAT (Ml  ) 

C  CONVERT  AMPERES  TO  KVA 
DO  60  1=1 f DAYS 
DO  50  J=2f25 

D ( I  • J )  " I F I X ( F L OA T ( D  < I y J ) ) * F A C TOR ) 

50  CONTINUE 
60  CONTINUE 
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KVADR 


C  CALCULATE  STANDARD  DEVIATION 
DO  90  J =2  s'  25 
S1=0 
S2®0 
lSi®o 
IS2=0 

DO  90  I®1 t DAYS 

IF  ( D  ( I » 1 )  . LT 4  0)  GOTO  70 
Sl=Sl+<D<IyJ>-M(lyJ> >**2 
IS1=IS1+ 1 
GOTO  80 

70  S2®82+ <  D  < I y  J  > -M ( 2  y  J  > ) *#2 

IS2=IS2+ 1 
80  CONTINUE 

S< 1 » J)*IFIX<SQRT< FLOAT < SI) /FLOAT < IS1-1) ) ) 

S ( 2  * J ) - IFIX  < SORT < FLOAT ( S2 ) /FLOAT < IS 2-1 ) ) > 

90  CONTINUE 

S ( 2  ? 1 ) 1 

C  WRITE  OUT  DATA  INTO  DATA  BASE  KVA. DBS 

.100  FORMAT  (  '  TYPE  IN  NAME  OF  FILE  WHICH  CONTAINS  DATA') 

110  FORMAT <A5> 

120  FORMAT <2A5?I5>  1.1  OX) 

130  FORMAT <2515) 

WRITE<7»120>  ( NAME2  < I ) y 1  =  1 y  2  > ? DAYS 
DO  140  I - 1 » 2 

WRITE  ( 7 » 130 )  <M<  I  >  J)  y  J~.l  ?25) 

WRITE ( 7  y 130 )  (S( ly J> y J-ly25) 

.1  40  CONTINUE 

DO  150  1=1 y DAYS 

WRITE ( 7  y 130)  (O(IyJ) ?J=1»25) 

1 5 0  CO N T I N U E 

CLOSE  ■'  UNIT=6  y  ACCESS®  '  SEQIN '  y  FILENAME  y  DEVICE*  ' DSK  '  ) 

CLOSE ( UNIT  =  7 r ACCESS® ' APPEND ' y  FILE® ' KVA . DBS ' t DEVICE® ' DSK  '  ) 
END 
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